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ABSTRACT 



Context. The gravitational lensing effect provides various ways to study the mass environment of galaxies. 

Aims. We investigate how galaxy-galaxy(-galaxy) lensing can be used to test models of galaxy formation and evolution. 

Methods. We consider two semi-analytic galaxy formation models based on the Millennium Run N-body simulation: the Durham 

model by Bower et al. (2006) and the Garching model by Guo et al. (2011). We generate mock lensing observations for the two 

models, and then employ Fast Fourier Transform methods to compute second- and third-order aperture statistics in the simulated 

fields for various galaxy samples. 

Results. We find that both models predict qualitatively similar aperture signals, but there are large quantitative differences. The 
Durham model predicts larger amplitudes in general. In both models, red galaxies exhibit stronger aperture signals than blue galaxies. 
Using these aperture measurements and assuming a linear deterministic bias model, we measure relative bias ratios of red and blue 
galaxy samples. We find that a linear deterministic bias is insufficient to describe the relative clustering of model galaxies below ten 
arcmin angular scales. Dividing galaxies into luminosity bins, the aperture signals decrease with decreasing luminosity for brighter 
galaxies, but increase again for fainter galaxies. This increase is likely an artifact due to too many faint satellite galaxies in massive 
group and cluster halos predicted by the models. 

Conclusions. Our study shows that galaxy-galaxy(-galaxy) lensing is a sensitive probe of galaxy evolution. 
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1. Introduction 

Gravitational lensing effects provide versatile tools for probing 
the matter distribution in the Universe. Galaxy-galaxy lensing 
(GGL), for example, is a statistical approach using lensing to 
obtain infor mation on the mass associated wit h individual galax- 
ies (see, e.g, iBartelmann & Schneider! l200ll) . This is achieved 
by dividing the galaxy population into lenses (foreground) and 
sources (background). The images of the sources are sheared 
due to the gravitational field of the foreground lenses and their 
surrounding mass. The image shearing is usually too small to 
be detected for individual source-lens galaxy pairs. Instead, the 
lensing effect is measured as a correlation between the observed 
image ellipticities and the lens positions. The signal obtained 
from averaging over many source-lens pairs can then be related 

to the average mass profiles o f the lenses. 

Since its first detection dBrainerd et al.|[T996h . GGL has 
been measured in many lar ge galaxy surveys (e.g.. Hoekstra 



et al.l2002t[Kleinheinrich et al.ll2006l:lMandelbaum et alJl2006t 
Simon et al. 2008; van Uitert et al-fcOlU and references therein). 
Schneider & Wattsl (I2005I) advanced GGL to galaxy-galaxy- 
galaxy lensing (G3L) by introducing third-order correlation 
functions that involve either configurations with two background 
sources and one lens galaxy (G+), or with two lenses and one 
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background source {Q). The latter measures the lensing signal 
around pairs of lens galaxies in excess of what one obtains by 
simply adding the average signals of two individual galaxies, 
and thus provides a measu re of the excess m atter profile about 
clustered lens galaxy pairs (ISimon et al l2012h . This G3L signal 
has been measured in th e Red sequence Clu ster Survey (RCS, 
iGladders & Yed2005l) bv lSimon et af] J20Q8h . who indeed found 
an excess mass about lens pairs with projected separation of 
250/r 1 kpc. The GG(G)L correlations can be converted to aper- 
ture statistics (which we utilize in this work), providing a conve- 
nient probe of the galaxy-matter power(bi-)spectra at particular 
scales. 



The galaxy-mass correlation as seen by weak lensing can 
also be studied theoretically by combining dark matter simu- 
lations with semi-anal ytic models (SAM) of galaxy evolution 
(IWhite & Frenklll99U iKauffmann etailll999l : ISpringel et all 
12001 ). In this approach, the dark matter halos of an A^-body sim- 
ulation of cosmic structure formation are populated with galax- 
ies. The properties of the galaxies are calculated by combining 
information on the halo merger trees of the underlying dark mat- 
ter simulation with an analytic model of the gas physics in galax- 
ies. The physical processes considered include gas cooling, star 
formation, metal enrichment, and feedback due t o supernovae 
and a ctive galactic nuclei. Using ray-tracing (e.g. iHilbert et al.l 
l2009l) . one can then simulate lensing observations of the result- 
ing galaxy distribution. 
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This paper provides a study of the second- and third-order 
galaxy-mass correlations in semi-analytic galaxy formation 
models as probed by lensin g via aperture statistics (ISchneiderl 
ll996USchneider et al.|ll998h. We consider two models based on 
the Millennium Run ([Springel et al.ll2.Q05b : the Dur ham model 
by Bower et al.l (12006 ) and the Garching model by iGuo et ail 
(1201 lh . We find that the predicted second- and third-order lens- 
ing signals differ between galaxies of different color and mag- 
nitude, but also between the different galaxy models. The dif- 
ferences between the models can be traced back to, among other 
things, different treatments of the satellite galaxy evolution. This 
illustrates that galaxy-galaxy(-galaxy) lensing can be a sensitive 
probe of galaxy evolution. 

The outline of the paper is as follows: Sect.|2]provides a brief 
account of gravitational lensing, aperture statistics, and their re- 
lation to correlation functions. Our lensing simulations and the 
method we use to measure aperture statistics (a fast method 
based on Fast Fourier Transforms) are described in Sect. |3] The 
results of these measurements for different subsets of galaxies, 
defined by redshift, luminosity or color, are presented in Sect. [4] 
The main part of the paper concludes with a summary and dis- 
cussion in Sect. [5] In the appendix, we briefly discuss shot-noise 
corrections for the aperture statistics. 



2. Theory 

2.1. Gravitational lensing basics 

The matter density inhomogeneities can be quantified by the di- 
mensionless density contrast 



Pm(x,x) ~ Pm(x) 
PmiX) 



(1) 



where p m (x,x) is the spatial matter density at comoving trans- 
verse position x and comoving radial distance x, and p m (x ) de- 
notes the mean density at that distance. To lowest order, the con- 
vergence K m for sources at comoving distance x is related to the 
matter densit y contrast 6 m by the projection along the line-of- 
sight by ("e.g. [Schneider et al.ll2006h 

3H 2 n m r* ,f K (x-x')fK(x')s m (f K (x')^x') 



hix) 



a(x') 



, (2) 



where K m describes the dimensionless projected matter density, 
Hq denotes the Hubble constant, Q. m the mean matter density pa- 
rameter, c the speed of light, /k(x) me comoving angular diam- 
eter distance, and a(x) = 1/(1 + z(x)) the scale factor at redshift 
Z(X). 

For a distribution of sources with probability density p$(x), 
the effective convergence is given by 

Km(&) = J dxPs(x) K m(&,X) 

dXg(x)dm(fK(x)&,X) w i th 



g(x) = 



2c 2 a(x) 



, , Jk(x' -x) 



(3) 



(4) 



Similar to the definition of the dimensionless matter density 
contrast S m , one can define the number density contrast 8 g of the 
lens galaxies as 



5g{x, X ) = 



p g (x,x) - p g (x) 
Pg(x) 



(5) 



where p g (x,x) is the number density of the lens galaxies at co- 
moving transverse position x and distance x> and p g is mean 
number density of lens galaxies at comoving distance x- Using 
the projected number density 



J dxfl(x)Pz(fK(x)#,X) 
and the mean projected number density 



(6) 



(7) 



the projected number density contrast for lens galaxies with dis- 
tance distribution p g (x) = n~ /lOf) Pg(x) can t> e computed by 



Kg(#) 



n g {d) - n g 



fd XPi 



(x)s»(f K (x)#,xl 



(8) 



2.2. Aperture Statistics 



Aperture statistics was originally introduced as a way to quantify 
the surface mass density tha t is unaffected by the mass sheet 
degeneracy (ISchneiderlll996h . The aperture mass is defined as a 
convolution, 



M ap (#;0) 



(9) 



of the convergence k and an axi-symmetric filter Ug{\&\) whose 
size is given by the scale 0, and that is compensated, i.e. 



d##t/ e (i?) = 0. 



(10) 



In this work, we use the filter introduced by ICrittenden et al.l 
(I2002h . 
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u(x) = — 
2n 



'4 



exp 
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(ID 
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Its (2-D) Fourier transform has a simple analytical form 

q2 £l 



Ug(€) 



U e (m)e™="—e-± B2e2 



u{k) = y ex P| — ]■ 



u{6€) with (13) 
(14) 



The filter falls off exponentially for § » 8. This makes the sup- 
port of the filter finite in practice. 

In analogy to the aperture mass M ap , one can define the aper- 
ture number count 



N(&- 8) 



U b (\#-&'\)kM'). 



(15) 



The aperture number count dispersion is related to the angular 
two-point correlation function 



g(02-0ll) = (*g(0l)*g(02)), 



(16) 
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and its Fourier transform, the angular power spectrum P gg (() of 
the lens galaxies, through 



(N 2 ){0) = (N(#; 6)N(#; 0)) 

= j d 2 #i U e (m) j d 2 § 2 U e (\& 2 \)w gg (\§ 2 -§ x \) 



(17) 



-—~u\e£)p gg {t). 

27T 



The function u(C0) features a sharp peak at 10 = V2. Thus, 
{N~)(0) provides a measurement of the corresponding power 
spectrum P gg (£) at wave numbers I ~ 1/0. 

Within the halo mo del framework of cosmic structure (e.g. 
ICoorav & Shefhll2002l) . (N 2 )(0) on small scales probes the 
distribution of the lens galaxies within individual dark matter 
halos. On large scales, (TV 2 ) provides a probe of the clustering 
of the host halos of the lens galaxies. 

Correlating M ap (6>) with N(0) yields 



(//M a p) (0) = (N(d; 0)M dp (#; 0)) 

= J" d 2 *i u e (m) J" d 2 § 2 u e (\# 2 \)w gm (\§ 2 ■ 



(18) 



me , 



with w gm (|#2 - = (%(^i)fm(^2)), whose Fourier transform 
is the cross-power spectrum of galaxies and convergence P gm . 
The galaxy-galaxy lensing aperture statistics (A/M ap ) probes the 

average matter profiles around lens ga laxies. 

A third-order aperture correlator (S chneider & Wattsll2005b 
is obtained by 

(/V 2 M ap ) (0) = (N(#; e)N(&; 0)M dp (#; 0)) 

= J" d^Ugim) J" d 2 # 2 U e (\# 2 \) J" d 2 # 3 U e (m) 
x(%(*ikg(«2>ft)) (19) 

= f tt^ f ^"(.Wiwmwm+hb 

J (2ttY J (In) 1 

X B ggm (t\,l 2 , -t\ —(2) 

where the last li ne contains the angular bi spectrum of the pro- 
jected quantities dSchneider & Wattsll2005h . 

<*i(fi)«a(f2)fc(* 3 )> = (2tt) 2 6 d (( 1 +e 2 + hW m (tuh,h). (20) 

On small scales, (N 2 M dp ) can teach us about the average mass 
distribution of halos hosting two lens galaxies. On larger scales, 
(AT 2 M ap ) also provides information on the higher-order cluster- 
ing of the host halos. 

2.3. Relative galaxy bias 

Clus ters and galaxies are biased tracer s of the matter d istribu- 
tion dKaisedll984lBardeen et al1ll986l : lMo et alJll996h . In the 
simplest conceivable non-trivial bias model, the bias can be ex- 
pressed as a linear deterministic relation between the galaxy den- 
sity contrast and the matter density contrast, 



S s {x,x) = b gm 5 m (x,x), 



(21) 



A more realistic assumption is that the galaxy bias is stochas- 
tic and depends on the time and spatial scale. At the two-point 
level, the bias may then be quantified by a scale-dependent bias 
factor and correlation factor. A description of the relation be- 
tween the galaxy and matter densities at the three-point level or 
higher requires additional, higher-order bias and correlation pa- 
rameters. 

Aperture stati s tics can be used to constrain the galaxy bias 
(ISchneiderl 119981: Ivan Waerbekd [19981: iHoekstra et al l 120021: 
ISchneider & Wattsll2005HSimon et alJl2007t Uullo et al.ll2012h 
For example, assuming a linear deterministic bias, (N 2 ) °c b 2 , 
(,/VM a p) oc b, and (N 2 M RV ) oc b 2 . For two lens galaxy samples 
with identical redshift distributions, but different bias parameters 
b\ and b 2 , one can then determine the relative bias b\jb 2 from 
the aperture statistics (N 2 / 2 ), {Ni/ 2 M ap ), and (N 2 / 2 M ap ) of the 
two lens samples by using any of 



bi l iNpW = (MM ap )(0) = (M 2 M. dp )(0) 
b 2 \<^ 2 2 )(0) {N 2 M. dV )(0) \{NlM dV ){9) 



(22) 



If the measured ratios of the aperture statistics change with 
scale 0, one can extend the idea to a measurement of a scale- 
dependent bias. For lens galaxy samples with narrow redshift 
distributions and deterministic bias, the above ratios still agree 
(roughly) when compared on the same scales. If t he galaxy bias 
is stochastic or non-linear (iDekel & Lahavl[l999l however, the 
ratios from the different statistics disagree even if measured on 
the same scale. In that case, the second- and third-order aperture 
statistics each contain valuable independent information on the 
secon d- and third-order bias of the galaxies (ISchneider & Watts! 
2005). 



3. Methods 

3.1. Lensing simulations 

For our analysis we u se the data obtained by ray-tracing through 
the Millennium Run dSpringel et al.l200"5l) . The Millennium Run 
(MR) is a large N-body simulation of structure formation in a flat 
ACDM universe with matter density Q m = 0.25, baryon density 
Ob = 0.045, dark-energy density Qa = 0.75, a Hubble con- 
stant Hq = MOOkm s _1 Mpc _1 with h = 0.73, and with a power 
spectrum normalization cr 8 = 0.9. It follows the evolution of 
N p ~ 10 10 dark matter panicles with mass m p = 8.6 x 10 8 /r' M Q 
in a cubic region of comoving side length 500/T 1 Mpc from red- 
shift z — 127 to the present. 

The simulation volume of the MR is large enough to in- 
clude massive rare objects, yet with sufficiently high spatial and 
mass resolution to resolve dark matter halos of galaxies. This 
allows the construction of merger trees of dark matter halos 
and subhalos within them. These merger trees have been used 
in various semi-analytic galaxy formation models to calculate 
the properties of galax ies in the sim u lation . Here we consider 
the Durh am model b y Bower et al.l (l2006h and the Garching 
model by iGuo et aTl (1201 1 ) 3 Both models have similar treat- 
ments of, e.g., gas cooling and star formation, but differ in var- 
ious details (see the original papers for a full description). The 
models have been adjusted to be consistent with a large num- 
ber of observations, in particular the luminosities, stellar masses, 



with a bias factor b gm that does not depend on time or scale, but 
only on the galaxy sample in question. 



1 We also considered an ea rlier incarnation of the Garching model 
bv lDe Luci a & Blaizot] te007l) . However, the differences between the 
results from two Garching models are minor. Th us, we concentrate th e 
discussion on the models bv lBower et ail ll2006h and lGuo et all d201lh . 
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morphologies, gas contents and correlations of galaxies at low 
redshift, but they have not been tuned to match galaxy proper- 
ties at higher redshift. We make use of the public Millenni um 
Simulations Databas^l dLemson & Virgo Consortiumll2006|) to 
obtain the properties of the galaxies predicted by the two mod- 
els. 

We em ploy the multiple-len s-plane ray-tracing algorithm de- 
scribed in Hilb ert et al ] (l2009h to calculate the light propaga- 
tion through the matter in the MR. We generate 64 simulated 
4x4 deg" fields of view. For each field, we calculate the con- 
vergence to sources at a number of redshifts on a regular mesh 
of 4096 2 pixels, as well as the apparent sky positions, redshifts, 
and magnitudes of the model galaxies from the Garching and 
Durham models. The galaxy properties are then used to select 
various subsamples of the full mock galaxy catalogs as lens pop- 
ulations for the GG(G)L. 



3.2. Computing aperture statistics 

We introduce a fast method to perform aperture statistics mea- 
surements on the 64 simulated fields. The statistic used here is 
built on two main components: (i) the pixelized convergence 
field K m (&) of the source galaxies on square meshes of 4096 2 
pixels, and (ii) the pixelized lens galaxy number density fields on 
meshes with the same geometry. The convergence fields are ob- 
tained directly from the ray-tracing algorithm. The galaxy den- 
sity fields are obtained by projecting the apparent position of the 
lens galaxies in the fields and counting the number of galaxies in 
each pixel. Finally, dividing by the mean number density of lens 
galaxies across all 64 fields results in the galaxy number density 
contrast /e g (#). 

We calculate the aperture statistics M ap (#; 9) and N(-&; 9) 
from K m and /c g on a grid by exploiting the convolution theo- 
rem, using Fast Fourier Transforms ( FFT, in particular the FFTW 
library by iFrigo & Johnsonl 120051) to carry out the convolu- 
tion in Eqs. (O and d 15b . To measure M ap (or N), we calcu- 
late the Fourier transforms of /c m (or /c g ) and UgQ&\). We then 
multiply the results in the Fourier space. Finally an inverse 
Fourier transformation gives M ap (or AO- The number of grid 
points in the field is finite (4096x4096 pixels). Therefore, a 
"Discrete Fourier Transform" is performed by using the Fast 
Fourier Transform (FFT) algorithm which reduces an 0(N 2 ) pro- 
cess to 0{N log 2 AO, with N being the number of points being 
Fourier transformed. Hence, FFT reduces the computation time 
immensely. 

Since M ap (#; 9) and A/($; 9) fields are not periodic, we ex- 
clude points closer than 4-9 to the field edges from the subsequent 
analysis. On the remaining points, we then calculate A/ 2 (#; 9), 
N(#; 0)M ap (#; 9), and N 2 {&; 6>)M ap (#; 9), and estimate (N 2 )(9), 
(NM ap )(9), and (A/ 2 M ap )(#) from these products by spatial av- 
eraging. We correct the estimates involving N for shot noise as 
described in AppendixlAl 

The G3L statistics can also be calculated from the shear field 
correspon ding to K m and the position s of the lens galaxies. In 
particular, ISchneider & Wattsl (120051) and ISimon et al"] (120081) 
showed that (A/ 2 M ap ) can be obtained as an integral over a three- 
point correlation function (3PCF). In order to check our proce- 
dure, we use 32 randomly selected simulated fields and calculate 
(A/ 2 M ap ) also with the latter method, by first calculating this 
3PCF with the help of a tree code. We note that, while this tree 
method is more flexible than the simple FFT-based method de- 
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Fig. 1. Upper panel: Aperture statistics (N 2 M. dp )(9) as a func- 
tion of filter scale 9 measured in the Garching model. The FFT 
method (squares) and the tree method (triangles) are compared 
for lenses at redshift z = 0.17 with m r < 22.5, stellar masses 
M* > 10 V M and convergence field of sources at redshift z = 
0.99. Error bars indicate the standard deviation of (A/ 2 M ap )(6>) 
for aperture radius 9 estimated across 32 fields. Lower panel: 
Average difference signal between the FFT method and the tree 
method. Again the error bars show the standard deviation of the 
mean (field variance of difference signal divided by "V31). 
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scribed above, in particular with regard to field boundaries and 
gaps, it is also considerably slower. 

For some of the individual fields, we find fairly large discrep- 
ancies between the results from the FFT and the tree method - in 
particular for fields with a large matter overdensity near the field 
boundaries. We can attribute these discrepancies to the differ- 
ent ways in which the three-point information is weighted in the 
two approaches. For example, a triplet of points near the bound- 
ary of the field enters the statistics in the tree method with the 
same weight as a similar triplet near the field center. In contrast, 
the FFT method, by excluding the stripe at the field boundary, 
assigns zero weight to such a triple. Hence, the results on indi- 
vidual fields can be quite different. 

Both methods are consistent, however, when averaging the 
results over many fields. Randomly selecting 32 simulated fields, 
we measure {N 2 M zp }(9) using the FFT method and the tree 
method. In the upper panel of Fig. [1] the outcomes of the two 
methods are compared, showing good agreement between the re- 
sults. The error bars, indicating the statistical error on the signal, 
tend to be smaller for the tree method than for the FFT method 
(for apertures larger than 2 arcmin), since the tree method makes 
better use of the fields' area. For example, for apertures larger 
than 20 arcmin, more than half of the field is not included in 
the FFT measurement. Consequently, the difference in scatter 
becomes more prominent on larger scales. The lower panel in 
Fig. Q] shows the field-by-field difference signal averaged over 
all fields. The difference between the methods is consistent with 
zero for 9 > 1 arcmin, but deviates from zero for 9 < 1 arcmin. 
This is due to a systematic underestimation of the signal in the 
tree method on small scales (ISimon et al . 2001. 
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Fig. 2. Redshift distribution of galaxies in the main lens samples 
(i.e. galaxies with redshifts 0.14 < z < 0.62, observer-frame r- 
band apparent magnitude m r < 22.5, and stellar masses M* > 



10 9 h M ) in the Garching and Durham model. 



4. Results 

4.1. Main lens samples 

In this section, we present results for the second- and third-order 
aperture cross-correlations and aperture number count disper- 
sion for the Durham and the Garching model in the 64 sim- 
ulated fields created from the Millennium Run. For simplicity, 
the background population is chosen to be located at z = 0.99. 
Unless stated otherwise, lens galaxies are selected to have red- 
shifts 0.14 < z < 0.62, observer-frame r-band apparent magni- 
tude m r < 22.5, and stellar masses M* > 10 9 /z 1 M . This yields 
8.5 x 10 6 lens galaxies in the Durham model and 8.7 x 10 6 galax- 
ies in the Garching model. The resulting redshift distributions for 
the lens populations are shown in Fig. [2] 

For this sample of galaxies, the aperture number count dis- 
persion (N 2 ) as a function of aperture radius is shown in the 
top panel of Fig. [3] The galaxy models clearly differ in the pre- 
dicted dispersion: the Durham model predicts an up to two times 
larger amplitude than the Garching model. A similar difference 
ha s been observed for the angular galaxy correlation function 
bv lKim et alJ d2009h , who attribute the discrepancy to too many 
bright satellites in the Durham model. However, as will be dis- 
cussed below, the Garching model also appears to suffer from 
problems with the modeling of the satellite population. 

The predictions for (JVM ap ), shown in the middle panel of 
Fig. [3] exhibit fairly large differences between the models, too. 
The higher values of (WM ap ) in the Durham model, especially 
for smaller angular scales, imply more massive lens halos on av- 
erage compared to the Garching model. The larger halo masses 
may also explain the higher clustering amplitude seen in (N 2 ). 
More massive halos host larger concentrations of galaxies and 
are themselves more clustered, which increases the correlation 
of the hosted galaxies on small and large scales. Another con- 
sequence is a larger third-order signal (A/ ,2 M ap ), which is con- 
firmed by the bottom panel of Fig. [3] 



4.2. Color-selected samples 

For a further analysis, we divide the main lens galaxy samples 
into groups selected by color. From observations, the color dis- 
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Fig. 3. Aperture number count dispersion (top panel), (NM ap ) 
(middle panel) and (N 2 M. dp ) (bottom panel) measurements in 
the Garching model and Durham model. 



tribution of galaxies is well characterized by a bimodal function 
dStrateva et al.l200"fl) . At low redshifts, this can be approximated 
by the sum of two Gaussian functions representing red and blue 
subpopulations of galaxies on the red and blue side of the color 
distribution, respectively. 

Following observations (e.g. iMandelbaum et al. 2006)), we 
split the main lens galaxy samples at observer-frame color u - r 
= 2.2 to obtain subsamples of red and blue lens galaxies. For 
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Fig. 4. Number of red and blue galaxies in the Garching model 
and the Durham model. Red galaxies are selected to have u — r> 
2.2 and blue galaxies are to have u - r < 2.2. 
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Fig. 7. Number of red and blue galaxies counted in the Garching 
model and the Durham model. Galaxies at each redshift are se- 
lected according to the color-cut in Fig. [6] 
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Fig. 6. The u — r color-cut at each redshift in the Garching and 
Durham models. 



the Durham model, we obtain 2.4xl0 6 red and 6.1xl0 6 blue 
galaxies compared to 4xl0 6 red and 4.7xl0 6 blue galaxies in 
the Garching model. The redshift distributions of the color sub- 
samples are shown in Figj4] The histograms show that the rela- 
tive numbers of red and blue galaxies in each redshift bin differ 
significantly between the models. 

The aperture statistics for these red and blue galaxy popu- 
lations in both models are shown in the left panels of Fig. 
In both models, red galaxies show higher signals than the blue 
galaxies. This trend is not surprising, since galaxies of different 
types follow different distribution patterns and clustering proper- 
ties. Red galaxies are expected to be found mainly in groups and 
clusters associated with strong clustering and large halo masses, 
whereas blue galaxies are mostly field galaxies with smaller ha- 
los and lower clustering. The plot also shows that galaxies in 
the two models show different clustering statistics in both red 
and blue populations. This may be a result of selecting different 
objects in the models. 

St udying the previous version of the Garch ing model based 
on the lDe Lucia & Blaizol (l2007h . iBetl (l2012h pointed out that 
the distributions of the observer-frame u - r colors are very dif- 
ferent in the Garching model and the Durham model. In partic- 



ular, galaxies appear redder in the Garching model than in the 
Durham model. We thus consider another way of dividing the 
galaxies into red and blue samples. We identify the minima of 
the bimodal distributions at each redshift. The positions of the 
minima are plotted in Fig. [6] clearly showing a large difference 
between the models in their color distribution. We then use the 
minima to separate red and blue galaxies. 

Fig. [7] shows the resulting redshift distributions of the red 
and blue subsamples. There are 4.2xl0 6 red and 4.3xl0 6 blue 
galaxies in the Durham model compared to 2.5xl0 6 red and 
6.2xl0 6 blue galaxies in the Garching model. Now the differ- 
ence between the models in the predicted numbers of red and 
blue galaxies is larger than for the case of a fixed color cut. 
This suggests that the redshift-dependent color cut selects very 
different objects in the two models. Surprisingly, the aperture 
statistics predicted by the two models are much more similar for 
the redshift-dependent color cuts than for the fixed color cut, as 
seen in the middle column of Fig . [5] The better agreement results 
from a decrease in the blue signals in the Durham model and an 
increase in the red signals in the Garching model. 

The agreement between the two models shown in the middle 
column of Fig. [5] indicates that although the redshift distribu- 
tions of red and blue galaxies differ, galaxies populate dark mat- 
ter halos in such a way to produce similar results. This agree- 
ment between the models is more prominently seen in (N 2 ) 
and (N 2 M ap ). Looking at (WM ap ), red Durham galaxies show 
a stronger signal on intermediate scales compared to the red 
Garching galaxies. This can happen, for example, if red galax- 
ies in the Durham model are mostly central galaxies populating 
large massive halos. On the other hand, this difference may also 
be a result of the distinct redshift distributions. As will discussed 
this difference is not seen when galaxies are restricted to come 
from a single redshift. 

Selecting lens galaxies at a single redshift amplifies the sig- 
nal for (AT 2 ) and (A/ ,2 M ap ) compared to a sample with a broad 
redshift distribution, where many of the projected galaxy pairs 
are at different redshifts and are therefore not correlated and 
suppress the overall signal. The third column of Fig. [5]displays 
the aperture measurements for lens galaxy populations selected 
from a single redshift slice around z = 0.17 with thickness 
Az = 0.02. Now the signals for red galaxies agree well between 
the two models. The agreement is not so good for blue model 
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Fig. 5. Aperture statistics for samples of red and blue galaxies in the Garching and Durham models. The left column shows the 
results for a fixed color-cut at u - r = 2.2 for galaxies with redshift between z = 0.14 and z = 0.62. The middle column displays the 
signals for galaxies between z — 0.14 and z = 0.62 separated using a redshift-dependent color cut. In the right column, galaxies are 
restricted to come from a single snapshot at redshift z = 0.17; accordingly error bars are larger. 
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galaxies, where the Garching model shows stronger signals on 
small and intermediate scales. In the halo-model language, blue 
Garching model galaxies at this redshift appear to live in more 
massive halos. 

Both the Durham and the Garching models predict a larger 
ratio between the clustering strength of red and blue galaxies 
than has been obtained in observations. A similar behavior was 
seen when consi dering the previous incarnati on of the Garching 
model based on iDe Lucia & Blaizol ( 20071). In particular, ou r 
results confirm the previous work of Ide la Torre et al. 
who compared the color-dependent projected two-point correla- 
tion function of a color sub sample of galaxies in t he VIMOS- 
VLT Deep Sur vey TV VPS: iLe Fevre et all 120051) and in the 
model based on IDe Lucia & Blaizotl d2007h 7 They showed that 
red galaxies in the semi-analytic models have stronger clustering 
amplitudes than the observed ones. They linked this discrepancy 
to an overproduction of bright red galaxies in the model. 

The different clustering strengths of red and blue galaxies 
show up very clearly in (N 2 ), (NM ap ), and (N 2 M ap ). The ratio 
of the clustering amplitude of the red and blue samples is related 
to their relative bias (Sect. l2~3V This ratio can be measured based 
on different aperture statistics measurements presented in Fig. |5] 
Assuming a simple linear deterministic bias, the relative bias and 
its uncertainty is calculated on aperture scales of 6 ~ 1 arcmin 
and 6 ~ 10 arcmin in the Garching and Durham models. The 
results are shown in Tab. Q] 



Table 1. The relative bias £> re d/£>biue based on different aperture 
statistics (Sect. 12.3b measured according to the right column of 
Fig. |5] on scales of 1 and 10 arcmin in the Garching and Durham 
models. The values are obtained assuming a linear deterministic 
bias model. 



(N 2 ) 



(NM ap ) (A/' 2 M ap ) 



1' 


Garching 


5.4±0.1 6.06±0.21 10.83±0.32 


Durham 


9.73±0.13 9.34±0.37 18.9+0.64 


10' 


Garching 


3.37±0.19 3.84±0.93 4.8±0.94 


Durham 


3.52±0.21 4.82+1.1 7.15±1.5 



The differences in the bias ratios measured from different 
statistics point out that a linear deterministic bias model is not 
sufficient to describe the relation between the galaxy and matter 
distribution on different scales. This relation may be described 
by scale-dependent stochastic bias. 

4.3. Magnitude-selected samples 

In this section, we present the measurements of the second- and 
third-order aperture statistics for lens galaxies in six different 
bins of r-band absolute observer frame magnitude M r . To elim- 
inate effects of possibly different redshift distributions on the 
signals, we restrict the redshift range of the lens galaxy popu- 
lation to one redshift slice at z = 0.17. The results for all magni- 
tude bins for the Garching (Durham) model are shown in the left 
(right) panels of Fig. [8] 

There are common trends seen for both models in the 
second- and third-order aperture statistics. For the brighter bins 
(-23 < M r < -20), the aperture signals decrease rapidly with in- 
creasing magnitude M r and filter scales 8. However, the Durham 
model predicts up to 200% higher (N 2 M- dp ), and (N 2 ) signals 
than the Garching model. Bright galaxies appear more clus- 
tered and on average to be located in more massive halos in the 
Durham model. 



For the fainter bins (-20 < M r < -18), the signals in- 
crease with decreasing luminosity. This increase is contrary to 
obse rvations of galaxy c lustering and GGL (see e.g. McBride 
et al. 1201 lh . where brighter galaxies show stronger clustering 
and larger lensing signals than fainter galaxies. 

In the Garching model, the faint magnitude bins are over- 
populated with satellite galaxies, many of which have no own 
subhalo (this occurs when a galaxy has been stripped of its own 
halo during a merger process with a larger halo). These galax- 
ies are abundant in massive halos, which contribute substantially 
to the (N M ap ) and (N 2 M. dp ) signal due to their large mass and 
stronger clustering. In the Durham model a similar trend is seen 
in (NM dp ) and (N 2 M dp ), indicating similar problems with the 
modeling of the satellite population in massive halos. 

The luminosity dependence of galaxy clustering has been 
studie d extensively with the aid of galaxy surveys. iLi et al.l 
(120071) compared the lumi nosity dependence o f the clustering 
of galaxies in the model of ICroton et al.l d2006l) to results from 
the Sloan Digital Sky Survey D ata Release Four (SDSS DR4; 
lAdelman-McCarfhv et al.ll2006l) . and found that the faint model 
galaxies show a stron ger clustering than SDSS galaxies. 

Kim et al. d2009h compared the galaxy clustering predicted 
by t he models of iBower et all (l2006l) , lDe Lucia & Blaizotl d2007h 
and iFont et al. I <l2008F to observed c lustering in the two- degree 
Field Galaxy Redshift Survey f2dF. IColless et al.ll200lh . They 
found that none of the models are able to match the observed 
clustering properties of galaxies in different luminosity bins. In 
particular, the Durham model shows a stronger signal than ex- 
pected, which could possibly be corrected, if the number of satel- 
lite gala xies in halos is re duced. 

Both ILi et al.l (120071) and iKim et alJ d2009h emphasize the 
problems of the galaxy models in predicting th e luminosity de- 
pendence of galaxy clustering. ILi et al.l d2007h showed that the 
number of faint satellite galaxies has to be reduced by 30 per 
cent (regardless of their h ost halo mass) to better match the 
observed galaxy clustering. IKim et al 1 (12001 showed that the 
fraction of s atellites declines with increasing luminosity (see, 
e.g. Fig 4 in IKim et alJl2009h in the host halo mass range of 
lO'V 1 M < M ha i < lO'V 1 M Q . Since the clustering strength 
depends strongly on the halo mass, satellite galaxies can affect 
the overall clustering amplitude. To investigate this they used a 
simple HOD model to show that satellite galaxies show a strong 
bias due to a strong two-halo clustering term. This indicates that 
satellite galaxies are preferentiall y found in massiv e halos which 
exhibits larger bias (see Fig. 5 in lKim et al.ll20Q9h . They argued 
that the results can be improved if the satellites are removed from 
massive halos by adding satellite-satellite merger processes in 
the models. 

Our results suggest that the Garching model shows a similar 
problem with faint satellite galaxies, though to a lesser degree. 
Indeed, we find that the amplitude of the aperture signals in the 
faint bins are completely dominated by satellite galaxies in both 
the Durham and the Garching model. Most of these faint satel- 
lites reside in massive group and cluster halos, which results in 
very large aperture signals on the scales considered in this work. 

5. Summary and discussion 

Through observations of galaxy-galaxy(-galaxy) lensing, valu- 
able information on the clustering properties of the galaxy 
and matter density field in the Universe can be obtained. 
Measurements of galaxy-galaxy lensing (GGL) can be used to 
infer information on the prop erties of dark matter h alos host- 
ing the lens galaxies (see, e.g JSchneider & Rixlll997| : Johnston 
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Fig. 8. Aperture measurements in the Garching model (left panel) and the Durham model (right panel) in 6 different r-band absolute 
magnitude, M r , bins. 
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et al. 120071: iMandelbaum et al.ll2008l) . Third-order galaxy lens- 
ing (G3L) can be used to infer information on the pr operties of a 
common da rk matter halo hosting two lens galaxies dSimon et alj 
I2008ll2012h . 

In this work, we study how the information from GGL and 
G3L aperture statistics can be used to test models of galaxy for- 
mation and evolution. We investigate two semi-analytic galaxy 
formation models based on t he Millennium Run N -body simu- 
lation of structure formation (ISpringel et al.ll2005l) : the Durham 
mod el bvlBower et al.l (120061). and the Garching model by_Guo 
et al. d201 lh . Using mock galaxy catalogs based on th ese models 
in conjunction with ray-tracing dHilbert et al.ll2009h . we create 
simulated fields of galaxy lensing surveys. From these simulated 
surveys, we compute the model predictions for the second- and 
third-order aperture statistics (N 2 ), (NM ap ), and (/V 2 M ap ) for 
various galaxy populations. 

We find that both semi-analytic models predict aperture sig- 
nals that are qualitatively similar, but there are large quantitative 
differences. The Durham model predicts larger amplitudes for 
most considered galaxy samples. This indicates that lens galax- 
ies in the Durham model tend to reside in more massive halos 
than lens galaxies in the Garching model. 

In both the Durham and the Garching model, red galaxies ex- 
hibit stronger aperture signals than blue galaxies, in qualitative 
agreement with observations. However, both models predict a 
larger ratio between the clustering strength of red and blue galax- 
ies than has been obta ined in observations. T hese findings cor- 
roborate the findings o flde la Torre et al. I dm ]]), who showed that 
red galaxies in the semi-analytic models have stronger cluster- 
ing amplitudes than red galaxies in observations. We argue that 
considering the amplitude ratio between the red and blue galax- 
ies and making comparison between the second- and third-order 
aperture statistics leads to the conclusion that the third-order 
bias differs from the second-order bias. In other words, third- 
order aperture statistics provides new information which cannot 
be obtained from the second-order statistics alone. The large am- 
plitude ratio between the clustering of red and blue galaxies in 
the models means a large relative bias of these galaxy popula- 
tions. Measuring the biasing of galaxies provides information 
on the relative distribution of galaxies and the underlying mat- 
ter distribution. We find that a linear deterministic bias model, 
even with scale-dependent bias parameters, is clearly ruled out 
by considering second and third-order aperture statistics for the 
simulated data. We expect that both statistics in combination will 
provide new information to constrain more advanced galaxy bi- 
asing models in the future. 

In addition to the different prediction for red and blue galax- 
ies, there are discrepancies between the predictions of the two 
models. For a fixed color cut at u — r = 2.2, the signals pre- 
dicted by the Durham model are larger than those predicted by 
the Garching model. If a redshift-dependent color cut is used in- 
stead, the prediction from the two models for the aperture signals 
become more similar. However, the models then strongly dis- 
agree about the total numbers and redshift distributions of blue 
and red galaxies. 

Both galaxy models predict that the aperture statistics de- 
crease with decreasing luminosity for brighter galaxies in accor- 
dance with observations. However, the models also predict that 
the signals increase again for fainter galaxies. This behavior is 
most likely an artifact related to too many faint satellite galax- 
ies in massive group and cluster halos predicted by the models. 
In fact, the fainter magnitude bins are completely dominated by 
satellite galaxies in both models. The problem appears more se- 



vere in the Durham model than in the Garching model, which 
differ in their treatment of satellite evolution. 

We plan to extent our treatment in future work to study how 
well galaxy bias models with scale-dependent stochastic bias 
can be constrained with second- and third-order galaxy lensing 
statistics. One important question is how much information can 
be obtained from G3L in addition to that obtained from GGL. 

Furthermore, we are looking forward to measurements of 
GGL and G3L signals in large ongoing and future surveys. The 
comparison of the observed signals and the signals predicted by 
galaxy models will help to identify shortcomings of the models 
and provide valuable hints for improvements in the models. This 
will also require a deeper understanding of the relation between 
the various details of the galaxy formation models and the pre- 
dicted galaxy lensing and clustering signals. 
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Appendix A: Shot-noise Correction 

Assume a realization (0. ), I = 1, . . . , N g , of a set of Ng galaxies with sky posi- 
tions fi. in a field M with area A and mean galaxy number density of h t = Ng/A. 
Assume that each galaxy position is distributed in the field according to an un- 
derlying 'true' number density field n s (&). The ensemble average of a quantity 
o({f) over all realizations reads: 



<0(0)> ■■ 



0(0). 



(A.l) 



For each realization, the random positions of the galaxies provide an esti- 
mate of the density field n g : 



This estimator is unbiased: 



(A.2) 



U N i J* 

1— f d 2 # r) n g (<' 



(A3) 



■■ nJff). 



Using a filter function U(&), we define a filtered density field N by: 
N(&, U)= \ d 2 &' U(& - &') n g (0') . (A.4) 



An estimator for the filtered field reads: 



N {r) (&; U)= f d 2 &' U(& - &')hf(&) = J^U(&- 



(A.5) 
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Its expectation value 



(N (r) (&; U)) r 



|-i£d 2 ^'n g (^)^-*f') 
f d 2 &' n g (&')U(& - &') 



(A.6) 



= JV(#; f/)- 

Consider the square of the filtered density 



N 2 (&;U) = f d 2 &' f d 2 &" U(&-&')U(&-&")n g (#')n g (#") 



(A.7) 



A naive estimator is provided by 



([JV w (#;{/)] 2 ) : 



x Y J u(&-# ( i r) )Y J u(&-# { j n ) 

I^I d2 ^l d2 ^ , " s( ^ ) " g < ) 

Y J ^-£d 2 ^ ) n,(& < p)U(&-4 



(A.8) 



iM)2 



A^g - 1) 



N 2 (&; U) + N(&; U 2 ). 



Hence, this estimator is biased. The first term of the last line is actually what is 
intended to be measured as aperture dispersion (up to a prefactor close to unity). 
The second term is due to shot noise arising from the Poisson sampling of the 
density field. An unbiased estimator of N~ is provided by 



Ng ■{W\&;U)] 2 -N (r) (&;U 2 )}- Ng 



N e - 1 



N„ ~ 1 



^ U(#-d l f ) )U(»-# ( -f > ). 



(A.9) 
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point- wise from the square of the former estimate to calculate the unbiased esti- 
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